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Abstract 



A stochastic control of the vibrational motion for a single trapped ion/atom 
is proposed. It is based on the possibility to continously monitor the motion 
Oi ■ through a light field meter. The output from the measurement process should 

be then used to modify the system's dynamics. 

> 

Pj ; 1 Introduction 

O 

\Q ' In recent years there has been an increasing interest on trapping phenomena and re- 

lated cooling techniques Ml. A few years ago, it has been shown that a single ion/atom 
can be trapped and cooled down near to its zero-point vibrational energy state 0|. 
On the other hand, the continuous measurement of quantum systems is particularly 

,G relevant at this time because experimental technology is now at the point where indi- 

vidual quantum systems can be monitored continuously 0] . With these developments 
it should be possible in the near future to control quantum systems in real time by 
using the results of the measurement in the process of continuous feedback. The 
possibility to control trapped particles, indeed, gave rise to new models in quantum 
computation M. In these systems the two lowest vibrational states are used for log- 
ical operations among quantum bits and the dominant source of decoherence is the 
heating of the vibrational motion B. Here, a way to control the position of a trapped 
ion/atom is studied. It results the possibility of a motional phase space uncertainty 
contraction corresponding to localise the particle even below the quantum limit. 



2 The Measurement Model 

Let us consider one vibrational mode for a trapped system together with a meter 
mode. The Hamiltonian would be of the following form 



1 email: mancini@mi.infn.it 



H = H m + fli/a f a + h x OX , (1) 

where H m is the free Hamiltonian of the meter, v being the oscillation frequency of the 
ion/atom in the trap, and a, a^ the lowering and rising operators for the vibrational 
states of the trap. In the last term of Eq.(|l]), x is the coupling constant, X = (a+af)/2 
is the dimensionless position operator, and O can be any meter operator obtained as 
linear combination of the -ff m 's eigenoperators O (but, for simplicity, in the following 
it is assumed as the sum of them). The practical implementation of such type of 
Hamiltonian will be discussed later on. 

The evolution equation for the whole density operator D is assumed to be 

t) = CD-i x [OX, D] + | (20-DO+ - O+Q-D - DO+Q-) , (2) 

where C describes the damped dynamics of the vibrational mode. Furthermore, the 
meter mode is considered heavily damped, so that the decay rate k is very large, (k 3> 
x), and the meter will almost always be in the lower eigenstate of H m . This allows 
the adiabatic elimination of the meter mode by means of a perturbative calculation 
in the small parameter x/ K i obtaining (see also Ref. g) 

D = [p <8> IT| + (-) [ipX®ILO- + h.c] + (-) [...] + ..., (3) 

where p — Tr TO D is the reduced density matrix for the vibrational motion, and II is 
the projector on the lower eigenstate of H m . 

Let us now soppose to be able to perform an homodyne-like measurement on the 
meter mode, i.e. to measure the observable O v — (0~e~ iv + + e lip )/2. Then, the 
current at readout apparatus will be []7) 

I(t) = 2 V K(O v (t)) c + JrjRZit) , (4) 

where r/ is the efficiency of the measurement process. The subscript c in Eq. (0) 
denotes the fact that the average is performed on the state conditioned on the results 
of the previous measurements and £(£) is a Gaussian white noise J7j. Therefore, due 
to the interaction in Eq.(l), one gets information on X by observing the quantity O v . 
From Eq.(pl), the relationship between the conditioned mean values results {O v ) c = 
(x/k)(X) c sin ip. 

However, the continuous monitoring of the meter mode modifies the time evolution 
of the whole system. In fact, the state conditioned on the result of measurement, 
described by the conditioned density matrix D c , evolves according to the following 
stochastic differential equation (considered in the Ito sense) 

D c = CD c -i X [OX,D c } + ^(20-D c O+ -0+0-D c -D c O+0-) 

+ ^m(e~ tV 0~D c + e^D c O+ -2(O v ) c D c ) . (5) 

If we adopt the perturbative expression (0) for D c in flq) and perform the trace over 
the meter mode, we get an equation for the reduced density matrix p c conditioned to 



the result of the measurement of the observable O v 
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£ Pc -^- [X, [X,/? C ]] + V?7X 2 MW (ie i(p p c X - ie-^Xp e + 2smtp(X) cPc ) . (6) 
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3 Applying the feedback loop 

Let us now consider the application of a feedback loop to control the dynamics of the 
vibrational mode. The continous feedback theory proposed by Wiseman and Milburn 
H is then invoked. One has to take part of the stochastic output current I(t), and 
feed it back to the vibrational dynamics in order to modify this. To be more specific, 
the presence of feedback modifies the evolution of the conditioned state p c (t)- It is 
reasonable to assume that the feedback effect can be described by an additional term 
in the master equation, linear in the photocurrent I(t), i.e. | 

\p&)] fb = y^r-icpc(t), (7) 

■qx 

where r is the time delay in the feedback loop and K, is a Liouville superoperator 
describing the way in which the feedback signal acts on the system of interest. 

The feedback term (|7|) has to be considered in the Stratonovich sense, since Eq. 
(m) is introduced as limit of a real process, then it should be transformed in the Ito 
sense and added to the evolution equation (|6|). A successive average over the white 
noise £(t) yields, in the limiting case of r — > 0, the master equation 

p = Cp - £ [ X , [X, p]} + K (ie^pX - ie-^Xp) + ^j^P- (8) 

The second term on the r.h.s. of Eq. (|S|) is the usual double-commutator term associ- 
ated to the measurement of X, it results from the elimination of the meter variables; 
the third term is the feedback term itself and the fourth term is a diffusion-like term, 
which is an unavoidable consequence of the noise introduced by the feedback itself. 

Due to the fact that the vibrational motion occurs at a frequency v of the order 
of MHz and the damping rate of the center-of-mass motion, 7, is usually very small 
|J , it seems reasonable [Q to use the generator L as that of quantum optical master 
equation at a nonzero temperature pd|| (indicating the number of thermal phonons 
with n = [exp(hv/kBT) — l] - *)- Moreover, since the Liouville superoperator K, can 
only be of Hamiltonian form S, it is taken as JCp = g \a — a\p\ /2, which means a 
driving term on the momentum, while g is the feedback gain related to the practical 
way of realizing the loop. Using the above expressions in Eq. (0) and rearranging the 
terms in an appropriate way, it is possible to get the following master equation: 

r L 

P = 77 ( N + !) (2a/9a f - a)ap - pa) a) + -N (2a) pa - aa) p - paa)) (9) 
M (2a ] pa) - a) 2 p - pa r2 ) M* (2apa - a 2 p - pa 2 ) - - sin^ [a 2 - a) 2 , p] , 
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Eq. (||) clearly shows the effects of the feedback loop on the vibrational mode a. The 
proposed feedback mechanism, indeed, not only introduces a driving term, but it also 
simulates the presence of a bath with nonstandard fluctuations, characterized by an 
effective damping constant T and by the coefficients M and N, which are given in 
terms of the feedback parameters M . An interesting aspect of this effective bath is 
that it is characterized by phase-sensitive fluctuations. 

Because of its linearity, the solution of Eq. (0) can be easily obtained in terms of 
the normally ordered characteristic function C(X, X*,t) [O. The stationary solution 
has the following form 



C(X, X* , oo) = exp 
where 



-C|A| 2 + ^(A*) 2 + ^*A 2 



(11) 



< = N + 9 ^ N ^ + ™*?}:^M 2 -, » = ^( C -N). (12) 

T z — g 2 sin ip gsmtp 

Under the stability conditions and in the long time limit (t — > oo) the variance of 
the generic quadrature operator Xg = {ae lB + a J( e~' ie )/2 becomes 

4(X e 2 ) = 1 + 2C + 2Re{ t ie 2t6 } . (13) 

For the position quadrature (9 = 0), it can be simply written as 

4(X 2 ) = 1 + 2n eff , 7i e// - C + Mm} • (14) 

In absence of feddback (g = 0) we have n e ff = n, otherwise n e ff can be smaller than 
n (the choice tp = —ir/2 turns out to be the best), providing a stochastic localisation 
in the position quadrature, i.e. a confinement. Depending on the external parameters, 
it can also be negative (but it is always n e ff > —1/2) accounting for the possibility 
of going beyond the standard quantum limit. This is a relevant result of the present 
feedback scheme since it is able to reduce not only the thermal fluctuations but even 
the quantum ones. 

4 The Practical Implementation 

Let us now consider a trapped ion/atom in a Lamb-Dicke regime, in the presence of 
a standing wave field, and let us study different ways to implement the above model 

El. 



A Hamiltonian resembling Eq. (|l|) can be obtained by considering an atom located 
at the node and in resonant condition for which [13] 

H = HAa z + hva f a + %x X , (15) 

with <j z , <j± the Pauli operators for two- level system, and A the-small-detuning be- 
tween the atomic and the standing wave frequency. The coupling constant is given 
by the Rabi frequency times the Lamb-Dicke parameter. Hence, in this case, the 
electronic degree of freedom plays the role of meter, i.e. a x = O and a± = ± . 
By exploiting the resonance fluorescence it could be possible to get the quantity 
S v = ((T-e~' lip + a+e lip ) /2 through homodyne detection of the field scattered by the 
ion along a certain direction |l4[ . In fact, the detected field may be written in terms 
of the dipole moment operator for the transition |— ) «-> |+) as E$ (t) — ^/rJK(T-(f) 
pjj |, where r\ is an overall quantum efficiency accounting for the detection efficiency 
and the fact that only a fraction of the fluorescent light is collected and superimposed 
with a mode-matched oscillator. 

Alternatively, the off resonant situation can be also exploited. The starting Hamil- 
tonian in this case implies the quantisation of the cavity field |ll| , and reads 

2e 2 
H = hAa z + hva'a + h-r-cx z Vbsin 2 (kX + cf>) (16) 

where b, w are the boson operators of the radiation mode, k is the dimensionless wave 
vector, e is the dipole coupling constant, and A the-large-detuning. 

By considering the internal atomic degree in the ground state and <p — 7r /4, the 
leading term of Eq. ( |16| ) is 

H = hva)a - h—b^b ( kX + 1/2) . (17) 

After linearization around the steady state of the cavity mode, the dynamics will be 
governed by an effective Hamiltonian formally identical to that of Eq. (|l|) 

H = h{-he 2 /A)b ] b + hva)a + h X YX , (18) 

where Y — (b+b')/2, and \ — —4/3ke 2 /A, with /3 the stationary value for the radiation 
amplitude (assumed real for semplicity). In this case the meter is represented by the 
cavity mode, i.e. O = Y, and + = a\ (D~ = a. The homodyne measurement of the 
light outgoing from the cavity allows to obtain the quantity Y v = (ae~ lip + a^e lip )/2 
as desired. 

Up to this point it has not considered how particular feedback Hamiltonians could 
be implemented. Of course it is important to be able to realize a term in the feedback 
Hamiltonian proportional to momentum. This is not so straightforward, but here 
possible ways in which it might be achieved are suggested. If the exact location of 
the trap is not an important consideration, then shifts in the position (being strictly 
equivalent to a linear momentum term in the Hamiltonian), are achieved simply by 



shifting all the position dependent terms in the Hamiltonian, in particular the trapping 
potential. This is a shift in the origin of the coordinates, and, being a virtual shift in 
the position, produces a term in the dynamical equation for the position proportional 
to the rate at which the trap is being shifted. When the experimental arrangement is 
such that the distance covered by the particle during the cooling is negligibly small 
compared to the trapping apparatus this may prove to be a very effective way of 
implementing a feedback Hamiltonian linear in momentum. Another method would 
be to apply a large impulse to the particle so that during one feedback time-step 
the particle is move the desired distance, and an equal and opposite impulse is then 
applied to reset the momentum. On the other hand, the use of laser pulses could be 
useful as well, since in accordance with the theory of laser cooling, the light exerts on 
the ion/atom a force proportional to its momentum. 

5 Conclusion 

Summarizing a stochastic control of the vibrational motion for a trapped particle via 
QND-mcdiatcd feedback has been proposed. In principle the model could be extended 
to the three dimensional case. The main limitation of such type of feedback is that 
it only acts on the measured (indirectly) variable, and it takes places only through a 
driving term in the variable conjugate to that measured. 

However, there are many ways in which the measurement signal may be fed back 
to affect the system. In general, at a given time, any integral of the measurement 
record up until that time may be used to alter the system Hamiltonian and affect 
the dynamics. Hence the above schemes could be improved. To this end a promising 
technique seems to be the feedback via estimation jig] . 

Acknowledgements The author acknowledges fruitful collaboration with P. Tombesi 
and D. Vitali at University of Camerino, Italy. 

References 

[1] see e.g., P. K. Ghosh, Ion Traps, (Clarendon, Oxford, 1995). 

[2] F. Diedrich, et ah, Phys. Rev. Lett. 62, 403 (1980). 

[3] H. Mabuchi, et ah, submitted to Appl. Phys. B, [eprint: |quant-ph/9805076[ . 



[4] J.I. Cirac and P. Zoller, Phys. Rev. Lett. 74, 4091 (1995). 

[5] D. DiVincenzo and B. Terhal, Phys. World 11, 53 (1998). 

[6] P. Tombesi and D. Vitali, Appl. Phys. B 60, S69 (1995); Phys. Rev. A 50, 4253 
(1994). 

[7] H.M. Wiseman, and G.J. Milburn, Phys. Rev. A 47, 642 (1993). 



[8] H.M. Wiseman and G.J. Milburn, Phys. Rev. Lett. 70, 548 (1993); Phys. Rev. 
A 49, 1350 (1994); H.M. Wiseman, Phys. Rev. A 49, 2133 (1994). 

[9] C. Monroe, et al., Science 272, 1131 (1996). 

[10] D.J. Wineland, et al, J. Res. Natl. Inst. Stand. Tcchnol. 103, 259 (1998). 

[11] see e.g., D. F. Walls and G. J. Milburn, Quantum Optics, (Springer, Berlin, 
1994). 

[12] S. Mancini, D. Vitali and P. Tombesi, submitted to Phys. Rev. Lett, [cprint 
|quant-ph/981002"2|. S. Mancini et al, in preparation. 



[13] C.A. Blockley, et al, Europhys. Lett. 17, 509 (1992). 

[14] J. T. Hoffges, et al, Opt. Comm. 133, 177 (1997). 

[15] K. Jacobs and A. Doherty, submitted [eprint [quant- ph/981200"4] 



